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Abstract 

The 3C-like proteinase of severe acute respiratory syndrome (SARS) coronavirus has been proposed to be a key target for structural- 
based drug design against SARS. Accurate determination of the dimer dissociation constant and the role of the N-finger (residues 1-7) 
will provide more insights into the enzyme catalytic mechanism of SARS 3CL proteinase. The dimer dissociation constant of the wild- 
type protein was determined to be 14.0 pM by analytical ultracentrifugation method. The N-finger fragment of the enzyme plays an 
important role in enzyme dimerization as shown in the crystal structure. Key residues in the N-finger have been studied by site-directed 
mutagenesis, enzyme assay, and analytical ultracentrifugation. A single mutation of M6A was found to be critical to maintain the dimer 
structure of the enzyme. The N-terminal octapeptide N8 and its mutants were also synthesized and tested for their potency as dimeriza¬ 
tion inhibitors. Peptide cleavage assay confirms that peptide N8 is a dimerization inhibitor with a K x of 2.20 mM. The comparison of the 
inhibitory activities of N8 and its mutants indicates that the hydrophobic interaction of Met-6 and the electrostatic interaction of Arg-4 
contribute most for inhibitor binding. This study describes the first example of inhibitors targeting the dimeric interface of SARS 3CL 
proteinase, providing a novel strategy for drug design against SARS and other coronaviruses. 

© 2005 Elsevier Inc. All rights reserved. 
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Human coronaviruses are major causes of upper respira¬ 
tory tract illness in humans. A novel form of coronavirus 
has been identified as the major cause of severe acute respi¬ 
ratory syndrome (SARS), a disease that had rapidly spread 
from southern China to several countries in 2003 [1,2]. Cor¬ 
onaviruses are members of positive-stranded RNA viruses 
featuring the largest viral RNA genomes known to date. 
The SARS coronavirus replicase gene encompasses two 
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overlapping translation products, polyprotein la 
(~450 kDa) and lab (~750 kDa), which are conserved 
both in length and amino acid sequence to other coronavi¬ 
rus replicase proteins. Polyprotein la and lab are cleaved 
by the internally encoded 3CL proteinase to release func¬ 
tional proteins necessary for virus replication. The SARS 
3CL proteinase is fully conserved among all released SARS 
coronavirus genome sequences and is highly homologous 
with other coronavirus 3CL proteinases. Due to the func¬ 
tional importance of SARS 3CL proteinase in the viral life 
cycle, it has been proposed to be a key target for structural- 
based drug design against SARS [3]. Virtual screening of 
chemical compound libraries has discovered potential 
inhibitors as described by various groups [4-7]. Several 
bifunctional boronic acid inhibitors were designed to target 
the serine cluster (Ser-139, Ser-144, and Ser-147) near the 
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catalytic residues [8]. Peptide analogs and organic com¬ 
pounds derived from known rhino virus 3C protease inhib¬ 
itors were also synthesized and some of them were found to 
be active towards SARS 3CL proteinase [9-14]. The crystal 
structures of the 3CL proteinases of several coronaviruses 
have been solved [3,15,16], all of which are dimeric. The 
enzyme was shown to follow a general base catalytic mech¬ 
anism [17]. As we have reported previously, the peptide 
cleavage assay shows that the specific activity for proteoly¬ 
sis decreases linearly with the decrease of enzyme concen¬ 
tration, suggesting that the dimer is the major form for 
biological activity and that the dimeric interface could be 
targeted for structural-based drug design against SARS 
3CL proteinase [18]. The extra helical domain of the SARS 
3CL proteinase was also shown to play an important role 
in controlling the association-dissociation equilibrium 
and regulating the activity and specificity of the enzyme 
[19]. However, there is no published report on the inhibitor 
design for the dimeric interface to date. 

The same strategy has been successfully used in the 
structure-based inhibitor design for human immunodefi¬ 
ciency virus 1 (HIV-1) protease and other viral enzymes 
[20-27]. HIV-1 protease is a homodimer of two identical 
99-residue subunits in which the active site is generated 
by the self-assembly of these subunits. The dimeric struc¬ 
ture is stabilized by an extended anti-parallel (3-sheet 
formed by the inter-digitations of the N-terminal (residues 
1^1) and C-terminal (residues 96-99) (3-strands of each 
monomer [22]. Peptides with the N- and C-terminal amino 
acid sequences have been designed and confirmed as dimer¬ 
ization inhibitors [21,27]. A large set of oligopeptides has 
been derived from the terminal fragments of HIV-1 prote¬ 
ase by computer modeling, and their inhibiting activities 
have been determined experimentally. The best inhibitor 
has an IC 50 value of less than 1 pM [20]. Cross-linked 
N- and C-terminal peptides also inhibit HIV-1 protease 
activity and reduce the fraction of dimers in solution as 
measured by size-exclusion chromatography [25]. Design 
of dimerization inhibitors has also been reported for many 
other viral enzymes, such as HIV reverse transcriptase, 
intergrase, herpes simplex virus ribonucleotide reductase, 
and DNA polymerase [24]. 

The crystal structure of SARS 3CL proteinase indicates 
that the N-finger fragment plays an important role in the 
dimerization and maintenance of the active form of the 
enzyme [16]. Peptides derived from this fragment may be 
used as dimerization inhibitors, as in the case of HIV-1 
protease and other viral enzymes. In the present paper, 
we have studied the key residues in the N-finger of SARS 
3CL proteinase by site-directed mutagenesis, activity test, 
and analytical ultracentrifugation. The N-terminal octa- 
peptide N8 and its nine mutants were also synthesized 
and tested for their inhibitory activities. The dimer dissoci¬ 
ation constants of the SARS 3CL proteinase and N-termi¬ 
nal mutants were determined by sedimentation velocity and 
equilibrium methods. This study provides important infor¬ 
mation for the structural-based inhibitor design targeting 


the dimeric interface of SARS 3CL proteinase, which pro¬ 
vides a novel strategy for drug design against SARS and 
other coronaviruses. 

Materials and methods 

The expression and purification of SARS 3CL proteinase. The C-ter¬ 
minal His-tagged 3CL proteinase was expressed and purified as described 
previously [18]. The non-His-tagged 3CL proteinase was expressed and 
purified as reported [17]. The R4E, K5A, and M6A mutants of SARS- 
CoV 3 CL proteinase were prepared with the QuikChange site-directed 
mutagenesis kit (Stratagene) using pET 3CLP-21h [17] as a template. The 
N-terminal deletion mutant, A7N, was constructed by inserting the PCR 
product which carried the Nhel/Hindlll site fragment with a deletion of 
the 6 codons coding for Gly-2 to Phe-7 into the Nhel/Hindlll-digested 
pET 21a DNA. The mutations were verified by nucleotide sequencing. The 
purification method used for each mutant was the same as that used for 
the wild-type proteinase and their purity was verified by SDS-PAGE to be 
greater than 95%. 

Peptide synthesis. The octapeptide inhibitor N8 and its mutants were 
synthesized by solid-phase peptide synthesis using the standard 9-fluore- 
nylmethoxycarbonyl/W/-butyl chemistry. Cleavage of the peptide from 
Rink resin and removal of all sidechain-protecting groups were achieved in 
trifluoroacetic acid solution. The crude peptide was purified by reversed- 
phase high-performance liquid chromatography (RP-HPLC, Elite P200II, 
Dalian, China) on a Zorbax Cl8 semi-preparative column (9.4 by 250 mm, 
Agilent) using water/acetonitrile gradients containing 0.1% trifluoroacetic 
acid. The purity of all peptides was measured by analytical RP-HPLC, and 
the identity was confirmed by matrix-assisted laser desorption/ionization 
time-of-flight mass spectroscopy. 

Colorimetric enzyme assay. The enzyme activity was measured by a 
colorimetric assay as reported before [17]. In short, 20 pi pNA substrate 
stock solution (2 mM Thr-Ser-Ala-Val-Leu-Gln-pNA water solution) was 
added to 180 pi 37 °C preheated reaction buffer (40 mM PBS, 1 mM 
EDTA, and 3 mM DTT, pH 7.3), which contained 2.8 pM enzyme. Col¬ 
orimetric measurements of enzymatic activity were performed in 96-well 
microtiter plates using a multiwell ultraviolet spectrometer (Spectra Max 
190, Molecular Device) at 390 nm. Each assay was repeated three times. 

Peptide cleavage assay. The proteolytic activity of the His-tagged 
SARS 3CL proteinase was determined using an HPLC-based peptide 
cleavage assay as previously reported [18]. The peptide substrates SOI: 
NH 2 -TSAVLQSGFRK-CONH 2 and SI2: NH 2 -SAVLQSGF-CONH 2 
were synthesized as described previously [18]. 

Zhang et al. [27] have established the Zhang-Poorman plot to distin¬ 
guish dimerization inhibitors from competitive inhibitors. The three major 
hypotheses in their model were: (i) the bioactive dimeric enzyme was in 
equilibrium with the inactive monomeric enzyme; (ii) the inhibitor only 
bound with the enzyme monomer; (iii) the substrate concentration in 
peptide cleavage assay was very low compared with K m of the dimeric 
enzyme. 

For the Zhang-Poorman plot, peptide S12 was used as substrate. 
Cleavage reaction mixture was incubated at room temperature and con¬ 
tained 400 pM substrate S12, 8.12-43.32 pM His-tagged SARS 3CL 
proteinase, no inhibitor or 800 pM peptide N8, and 57 pM DTT in 40 mM 
Tris-HCl buffer, pH 7.3, at room temperature. Reaction was stopped after 
5 or 10 min by adding 0.1% trifluoroacetic acid aqueous solution and 
analyzed as described previously. 

The inhibitor peptide was dissolved in DMSO and premixed with the 
His-tagged proteinase for 1-2 h. Cleavage reaction mixture (400 pM 
substrate SOI, 2.43 pM His-tagged SARS 3CL proteinase, 800 pM 
inhibitor peptide, and 57 pM DTT in 10% DMSO, 40 mM Tris-HCl 
buffer, pH 7.3) was incubated at room temperature. Aliquots of reac¬ 
tions were removed every 20 min within 80 min and the reaction was 
stopped by the addition of 0.1% trifluoroacetic acid aqueous solution, 
and the reaction mixture was analyzed by RP-HPLC as described above. 
k ca fK m was determined by plotting the substrate peak area using the 
equation below: 
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InPA = C - {k cat /K m ) app c E t, (1) 

where PA is the peak area of the substrate peptide, c E is the total concen¬ 
tration of His-tagged 3CL proteinase, and C is an experimental constant. 
The inhibitory activities of the peptides were estimated using the same 
dimerization inhibitor model as in the Zhang-Poorman plot [27]. The 
dimerization inhibition constant K x was thus calculated using the following 
equation: 

= 1+ !?’ (2) 

V (^cat/A m ) app J K x 

where ( k cat /K m ) app _o is the apparent k cat /K m without inhibitor and [7] is the 
concentration of the inhibitor. 

Analytical ultracentrifugation analysis. Sedimentation equilibrium and 
velocity experiments were conducted on a Beckman Optima XLA ana¬ 
lytical ultracentrifuge equipped with absorbance optics. An An60Ti rotor 
and standard six-sector equilibrium centerpieces were used. The freshly 
prepared wild-type and mutational SARS 3CL proteinase was further 
purified and buffer-exchanged using a gel filtration column, Superdex 75 
10/300 GL (Amersham Bioscience), into sedimentation buffer (40 mM 
phosphate buffer, 100 mM NaCl, and 0.5 mM EDTA, 0.5 mM DTT, pH 
7.3). The molar extinction coefficient at 280 nm (1.04 mg/cm 2 ), density of 
the sedimentation buffer (1.005 g/ml), partial specific volume (0.723 ml/g), 
and molecular weight of the monomer (33 914 Da) were calculated based 
on its amino acid composition using the program SEDNTERP (http:// 
www.bbri.org/rasmb/rasmb.html). For sedimentation velocity experi¬ 
ments, 380 pi samples (concentrations between 0.1 and 4 mg/ml) and 
400 pi reference solutions were loaded into cells. The rotor temperature 
was equilibrated at 20 °C and rotor speeds of 60,000 rpm. Absorbance 
scans at 230, 280 or 290 nm were collected at a time interval of 4 min. Data 
were analyzed with the software Sedfit version 8.9 g [28,29]. For the sed¬ 
imentation equilibrium experiments, 110 pi samples and 120 pi reference 
solutions were loaded into nitrogen-flushed cells, followed by degassing 
and a further nitrogen flush prior to sealing. The protein was equilibrated 
for data collection at 20 °C and three rotor speeds (15,000, 20,000, and 
25,000 rpm). Once equilibrium was reached (typically 24-32 h), absorption 
data were collected at 280 nm, using a radial step size of 0.001 cm, and 
recorded as the average of 10 measurements at each radial position. To 
determine the baseline values in the cell, at the end of the data collection 
time the rotor speed was increased to 42,000 rpm for 8 h, and the absor¬ 
bance of the depleted meniscus was measured. Dissociation constants were 
determined by fitting a monomer-dimer equilibrium model using the 
Origin-based data analysis software for Beckman XL-A/XL-I (Beckman 
Instruments, Beckman Coulter, Fullerton, CA). Data from different con¬ 
centrations and speeds were combined for global fitting. 

Results and discussion 

Dissociation constant of SARS 3CL proteinase dimer 

The sedimentation experiments study the aggregation 
state of the enzyme in the native solution condition, and 
no dilution effect exists as compared to gel filtration. 
AUC is considered more reliable than gel filtration in deter¬ 
mining protein aggregation states. Sedimentation velocity 
provides hydrodynamics information about the molecular 
size distribution and conformational changes, whereas sed¬ 
imentation equilibrium can perfectly characterize the ther¬ 
modynamics of the self-association systems and determine 
the dissociation constants [30,31]. Recently, Peter Schuck 
has developed a method for the analysis of protein self-as- 
sociation by sedimentation velocity experiments [29,32]. In 
order to get a reliable result, we performed both sedimen¬ 
tation velocity and equilibrium experiments to determine 


the dissociation constant of SARS 3CL proteinase. We 
determined the dissociation constant of the wild-type pro¬ 
teinase with no tags. 

Sedimentation velocity experiments on the wild-type 
enzyme were conducted at nine enzyme concentrations. 
The resulting c(s) distribution profiles show concentra¬ 
tion-dependent peaks at positions between the monomer 
and dimer (Fig. 1A). This clearly indicates the fast revers¬ 
ible association of the SARS 3CL proteinase. The isotherm 
analysis of s w (c) is shown in Fig. IB based on the integra¬ 
tion of the differential sedimentation coefficient c(s) distri¬ 
butions [29]. By fitting to the monomer-dimer 
equilibrium in the Sedphat program, the dissociation con¬ 
stant was determined to be 22.9 pM, with the monomer 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 


Sedimentation Coefficiency (S) 



Fig. 1. (A) SARS 3CF proteinase sedimentation coefficient distribution 
c(s) of nine SARS 3CF proteinase concentrations from 0.11 to 5.5 mg/ml. 
(B) Dimer dissociation constant analyzed by the isotherm analysis of 
weight-average sedimentation coefficient as a function of concentration. 
The raw data (circle) were obtained by the integration of the differential 
sedimentation coefficient c(s) distributions of nine enzyme concentrations. 
Sedphat program was used to fit the data to the monomer-dimer 
equilibrium (solid line). The resulted monomer sedimentation coefficient 
Si = 2.7 s, dimer sedimentation coefficient S 2 = 4.0, and the dissociation 
constant K d — 22.9 pM. 
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sedimentation coefficient of Si = 2.7 s, the dimer sedimen¬ 
tation coefficient of S 2 = 4.0 s (Fig. IB). 

The sedimentation equilibrium experiment was per¬ 
formed with three different protein concentrations (0.2, 
0.3, and 0.5 mg/ml) and three running speeds (15,000, 
20,000, and 25,000 rpm). Fig. 2 shows the result of a global 
fit of all the data to the monomer-dimer equilibrium with 
the Beckman Origin-based data analysis software. The dis¬ 
sociation constant for the monomer-dimer self-association 
model was determined to be K d = 14.0 pM, which is in 
agreement with the value from the sedimentation velocity 
experiment. Since the sedimentation velocity method is 
only a good initial exploration to characterize the self-asso- 
ciation of proteins, we take the final K d to be 14.0 pM from 
the global fit of sedimentation equilibrium data. 

There are several reports on the quantitative measure¬ 
ments of the dimer dissociation constant of the SARS 
3CL proteinase by analytical ultracentrifuge or enzyme 
activity-dependent assays. But the published dissociation 
constants have large variations from different groups 
and/or using proteins expressed from different constructs. 
Chou et al. [33] reported a dissociation constant of 
190 ± 14 nM by using the sedimentation velocity method 
with the C-terminal His-tagged form of SARS 3CL pro¬ 
teinase. Hsu et al. [34] reported a dissociation constant of 
0.35 nM by using the sedimentation velocity method on 
the wild-type protein with no tags. Our sedimentation equi¬ 
librium experiments give a dissociation constant of 
14.0 pM for the wild-type protein with no tag but one extra 
alanine residue at the N-terminus. Generally speaking, the 
sedimentation equilibrium test can provide first-principle 
analysis of the dissociation constants compared to the sed¬ 
imentation velocity method [30,31]. 



Fig. 2. Sedimentation equilibrium data for wild-type SARS 3CL protein¬ 
ase performed at 20 °C with three different rotor speeds (15,000, 20,000, 
and 25,000 rpm) and three different enzyme concentrations (0.2, 0.3, and 
0.5 mg/ml). The bottom panel is the nine equilibrium data and the best fit 
(solid line); the upper panel is the residuals’ plot. The dissociation constant 
( K d ) was determined to be 14.0 pM. 


Key residue analysis of the N-finger of SARS 3CL 
proteinase 

The crystal structure of SARS 3CL proteinase suggests 
that the N-finger fragment contributes significantly to the 
enzyme dimerization [16]. Our experiments showed that 
the N-finger deletion mutant A7N could not form any 
dimer (Fig. 3B), which agrees with the previous report that 
the N-finger deletion mutant is in monomeric form in solu¬ 
tion [35]. However, the seven residues in the N-finger con¬ 
tribute differently to the dimer stability. From Fig. 4, we 
can see that the side chains of Phe-3 and Lys-5 point away 
from the dimer interface to the interior of the same mono¬ 
mer. Ser-1, Gly-2, and Ala-7 do not have remarkable inter¬ 
actions. The two most important residues among the seven 
are Arg-4 and Met-6. The guanidyl group of Arg-4 in one 
monomer forms a salt bridge with the side chain of Glu-290 
in the other monomer, which is one of the major interac¬ 
tions between the two monomers. Furthermore, the side 
chain of Met-6 forms hydrophobic interactions with the 
side chains of the Phe-140 and Tye-126 in its partner. This 
suggests that residues Arg-4 and Met-6 are the two possible 
key residues for the N-finger role in dimer stability. How¬ 
ever, Chou et al. [33] reported that the mutation of R4A 
only partly affects the dimerization of SARS 3CL protein¬ 
ase and keeps most enzymatic activity, which means that 
Arg-4 might not be the dominating factor to the dimer 
interface. In this study, R4E and M6A mutants were shown 
to have very weak dimerization according to sedimentation 
experiments (Figs. 3C and E) and have no detectable 
enzyme activity either (Table 1). Substituting the positively 
charged Arg-4 with a negatively charged Glu makes it 
impossible to form a salt bridge with the Glu-290 in the 
other monomer and even cause charge repulsion. The 
mutation of M6A suggests that the hydrophobic interac¬ 
tions between the side chain of Met-6 and Phe-140, Tye- 
126 in the other monomer contribute significantly to the 
stability of the dimer. 

Our study confirms that the N-finger plays an important 
role in stabilizing the dimer of SARS 3CL proteinase. The 
N-finger of one monomer is important for the enzyme 
activity of the other monomer [15,16]. So it is easy to 
understand the phenomenon of no dimer, no activity. 
The mutation of R4A kept most of the enzymic activity 
and did not affect the dimerization severely [33], which 
means that the electrostatic interaction contributed by the 
salt bridge between Arg-4 and Glu-290 may not be vital 
to the dimerization. However, the alanine substitution of 
Met-6 would be lethal to the dimer. In fact, Arg-4 is not 
conserved and can be substituted by lysine or valine in sev¬ 
eral other coronavirus 3CL proteinases, while Met-6 is 
highly conserved with Leu substitution in very few cases 
[15]. As a single Met-6 to Ala mutation makes the enzyme 
to be mostly a monomer with no detectable activity, the 
hydrophobic effect of Met-6 should play a key role in the 
dimerization. Since dimerization of SARS 3CL proteinase 
is necessary for the enzyme activity, the interface of the 
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Fig. 3. Sedimentation velocity experiments of SARS 3CL proteinase (WT and mutants). Each enzyme solution was run with both high and low 
concentrations at 60,000 rpm and 20 °C. (A-E) Show the continuous sedimentation coefficient distribution c(s) of wild-type enzyme and mutants. (A) Wild 
type; (B) A7N; (C) R4E; (D) K5A; (E) M6A. 



Fig. 4. Electrostatic interaction between one monomer of SARS 3CL 
proteinase and the N-terminal octapeptide of its partner. Figure was 
drawn by Molmol using its simple charge parameters and the 1UJ1 
coordinates from PDB. Met-6 in the N-terminal octapeptide was 
accommodated in a hydrophobic pocket and Ser-1 and the sidechain of 
Arg-4 are placed in two negatively charged pockets. 

dimer may be used as a new target for structure-based drug 
design [18,19]. From the N-terminal mutational analysis, 
we propose that a hydrophobic group and an acidic moiety 


Table 1 


Activity of the mutational SARS-CoV 3CL proteinase 


Enzyme 

Relative catalytic activity 3 (%) 

WT 

100 

A7N 

<0.01 

R4E 

1 

K5A 

70 

M6A 

<1 


a The concentration of the enzyme is 2.8 pM. 


that can form strong hydrogen bonds with ARG-4 are 
essential to design strong dimer interface binding 
inhibitors. 

The N-finger fragment octapeptide N8 is a dimerization 
inhibitor of SARS 3CL proteinase 

The crystal structure of SARS 3CL proteinase indicates 
that the N-finger fragment contributes significantly to the 
enzyme dimerization [16], which may act as a dimerization 
inhibitor of the enzyme, as in the case of HIV proteinase 
and other viral enzymes. To determine whether the N-fin¬ 
ger fragment is inhibitory towards SARS 3CL proteinase, 
we have synthesized an octapeptide N8, corresponding to 
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the N-terminal residues 1-8 of the enzyme. The peptide 
inhibitor (and other peptide mutants) was premixed with 
the His-tagged enzyme before peptide cleavage assay for 
1-2 h. The peptide N8 has obvious inhibitory activity at 
millimolar level. 

As shown in Table 2, the changes in reaction rate were 
caused mainly by the changes in apparent /± at . The appar¬ 
ent K m of the His-tagged proteinase changes little at differ¬ 
ent concentrations of N8, indicating that the N8 peptide 
works as a non-competitive inhibitor. In other words, the 
N8 peptide binds to a site distinct from the substrate pock¬ 
et or to a different enzyme conformation. 

Zhang et al. [27] have established the Zhang-Poorman 
analysis to distinguish dimerization inhibitors from com¬ 
petitive inhibitors by plotting c E / ^/k exp vs. ^/k QXp in their 
research about the dimerization inhibitor of HIV protease. 
In the Zhang-Poorman analysis, a dimerization inhibitor 
alters the y-intercept value by a factor of [T\/K E but has 
no effect on the slope [20,27]. Dimerization inhibitors are 
easily identified using the Zhang-Poorman plot as they 
provide lines that are parallel with that from the protease 
alone, while competitive and non-competitive inhibitors 
give non-parallel lines [36]. Zhang-Poorman plot has been 
successfully applied in many protease dimerization inhibi¬ 
tor studies [20,23,25,27,36]. Fig. 5 shows that the line in 


Table 2 


The apparent K m and k cat of the SARS 
concentrations of inhibitor N8 

3CL proteinase 

at different 

Concentration of N8 (mM) 

0 

0.4 

0.8 

Apparent K m d (mM) 

1.00 ±0.03 

0.96 ± 0.05 

1.01 ±0.07 

Apparent k c at a (min -1 ) 

18.0 ±0.4 

12.7 ±0.6 

8.2 ±0.6 


a The concentration of the enzyme is 1.22 |aM. 



Fig. 5. Zhang-Poorman plot of C-terminal His-tagged 3CL proteinase. 
The lower line (■, without inhibitor) indicates that the enzyme is 
competitive with Zhang-Poorman’s model, in which only the dimer of the 
enzyme is biologically active. The upper line (O, with 0.8 mM peptide N8) 
is parallel with the lower line, indicating that the inhibitor N8 is a 
dimerization inhibitor of the proteinase, which prevents the enzyme 
proteolysis activity by disturbing the natural proteinase dimerization. 


the presence of 0.8 mM peptide N8 is parallel with that 
in the absence of inhibitor. This strongly suggests that 
the inhibitor peptide N8 is a dimerization inhibitor of 
SARS 3CL proteinase, which prevents the enzyme proteol¬ 
ysis activity by disturbing the dimerization of the 
proteinase. 

SARS 3CL proteinase has been proposed to be a key 
target for structural-based drug design against SARS, 
because of its functional importance in the viral life cycle. 
However, all published reports on inhibitor design to target 
the proteinase are based on the substrate binding pocket 
[4-7,9-14] and the serine cluster near it [8]. There is no 
report on the inhibitor design targeting the dimeric inter¬ 
face to date. The octapeptide N8 is the first inhibitor 
designed to bind to the dimeric interface of coronavirus 
3CL proteinase, thus providing a novel strategy for drug 
design against SARS and other coronaviruses. 

The different inhibitory activities of N8 mutants 

To further investigate the inhibition mechanism of N8, 
nine mutants were synthesized and their inhibitory activi¬ 
ties were determined (see Table 3). According to our muta¬ 
tional study, M6A and R4E mutants disrupted the dimer 
structure. For the N8 peptide inhibitor, mutation of Arg- 
4 to Ala (peptide N82) completely abolishes the inhibitory 
activity. Other single mutations of Arg-4 and Met-6 (pep¬ 
tide N81 and N84) also decrease N8’s inhibition activity 
significantly. These results correlate well with the mutation¬ 
al study. 

It is interesting that peptide N86 containing double 
mutations (Arg-4 to Ala and Met-6 to Ala) shows some 
inhibitory activity, even though it is predicted to be inactive 
based on the results from the single mutants. Since the N- 
finger fragment adopts an extended conformation in the 
crystal structure, the side chain of Lys-5 points outwards 
to the dimer interface and does not participate in the 
enzyme dimerization. The mutation of Lys-5 to Ala (pep¬ 
tide N83) only slightly decreases the inhibiting activity of 
the N-finger fragment, indicating that Lys-5 contributes lit¬ 
tle to the inhibitor binding. However, in peptide N86, the 
peptide backbone may turn around with the Lys-5 taking 
the place of Arg-4 and interacting with the enzyme, due 
to the loss of the conformation restriction offered originally 
by Arg-4 and Met-6. This is confirmed by the triple-muta¬ 
tion peptide N87, which shows little inhibition activity. 

In summary, we have carried out a mutational study on 
the N-finger of SARS 3CL proteinase and determined the 
dimer dissociation constants for the wild-type protein and 
the mutants using sedimentation velocity and equilibrium 
techniques. The Met-6 residue was found to be essential 
for maintaining the dimer structure, while Arg-4 contribut¬ 
ed less. The N-finger octapeptide N8 has been synthesized 
and confirmed to be a dimerization inhibitor with a K[ of 
2.20 mM. Several N8 peptide mutants were also synthe¬ 
sized and tested for their inhibition activity. The hydropho¬ 
bic interaction of Met-6 and the electrostatic interaction of 
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Table 3 


The inhibiting activities of N8 and its mutants 



Sequence 3 

(^'cat/X m )app (mM min ) 

Inhibitory (%) 

Experimental pKi 

No inhibitor 


23.6 ±0.9 

0 


N8 

sgfrkmaf-nh 2 

12.7 ±0.2 

46.2 

2.66 

N81 

sgfekmaf-nh 2 

17.2 ±0.8 

27.1 

2.33 

N82 

sgfakmaf-nh 2 

22.6 ±0.6 

4.2 

1.44 

N83 

sgframaf-nh 2 

14.0 ±0.3 

40.7 

2.57 

N84 

sgfrkaaf-nh 2 

19.4 ± 1.7 

17.8 

2.11 

N85 

SGFEAMAF-NHo 

18.0 ±0.8 

23.7 

2.26 

N86 

sgfakaaf-nh 2 

18.2 ±0.8 

22.9 

2.24 

N87 

SGFEAAAF-NH, 

21.3 ±0.2 

9.7 

1.82 

N88 

Ac-SGFRKMAF-NH? 

13.6 ±0.6 

42.4 

2.60 

N89 

Ac-SGFEAAAF-NH 2 

21.1 ± 1.0 

10.6 

1.86 


a The mutant residues are underlined. 
b The concentration of the enzyme is 2.43 pM. 


Arg-4 in N8 play important roles in inhibitor binding. Pep¬ 
tide N8 is the first inhibitor designed to target the dimeric 
interface of coronavirus 3CL proteinase. Our current study 
provides a novel strategy for drug design against SARS 
and other coronaviruses. 
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